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ABSTRACT 

Merging clusters of galaxies are unique in their power to directly probe and place limits on the 
self-interaction cross-section of dark matter. Detailed observations of several merging clusters have 
shown the intracluster gas to be displaced from the centroids of dark matter and galaxy density by 
ram pressure, while the latter components are spatially coincident, consistent with collisionless dark 
matter. This has been used to place upper limits on the dark matter particle self-inteaction cross- 
section of order 1 cm 2 g _1 . The cluster Abell 520 has been seen as a possible exception. We revisit 
A520 presenting new HST ACS mosaic images and a Magellan image set. We perform a detailed 
weak lensing analysis and show that the weak lensing mass measurements and morphologies of the 
core galaxy-filled structures are mostly in good agreement with previous works. There is however 
one significant difference - we do not detect the previously claimed "dark core" that contains excess 
mass with no significant galaxy overdensity at the location of the X-ray plasma. This peak has been 
suggested to be indicative of a large self-interaction cross-section for dark matter (at least ~ 5a larger 
than the upper limit of 0.7cm 2 g _1 determined by observations of the Bullet Cluster). We find no such 
indication and instead find that the mass distribution of A520, after subtraction of the X-ray plasma 
mass, is in good agreement with the luminosity distribution of the cluster galaxies. We conclude that 
A520 shows no evidence to contradict the collisionless dark matter scenario. 
Subject headings: Gravitational lensing - Galaxies: clusters: individual: A520 - dark matter 



1. INTRODUCTION 

Merging clusters of galaxies provide a unique oppor- 
tunity to study properties of dark matter. During a 
cluster merger, the cluster galaxies are effectively col- 
lisionless particles, slowed only by tidal interactions, 
while the ionized X-ray emitting plasma clouds are 
affected by ram pressure as they pass through each 
other. The ram pressure causes the plasma clouds 
to slow down, and shortly after each collision in the 
merger process the X-ray plasma clouds will be found 
betw een the major concentrations o f cluster galaxies 
(e.g. iRoettiger. Loken. fc Burns! Il997j ). Any dark mat- 
ter present would be located in the vicinity of the cluster 
galaxies, provided the dark matter does not a large self- 
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interaction cross-section. Because the X-ray plasma in 
a cluster makes up ~ 12% of the mass of a rich cluster 
(jAllen et al.l 120021: IVikhlinin etall I2006D while the stel- 
lar mass in the clus ter galaxies is less th en 1% of the 
mass of the cluster ([Kochanek et al.l 120031) , one has the 
situation that shortly after each collision the bulk of the 
baryonic matter is spatially displaced from the the bulk 
of the total mass, and thus from the largest gravitational 
potentials in the cluster. By comparing the positions 
of the X-ray plasma clouds to the centers of the grav- 
itational potential, as revealed using gravitational lens- 
ing measurements, one can place measurements on the 
properties of dark matter and test alternative theories of 
gravity, which one could use to replace some or all of the 
dark matter in galaxies and clusters. 
The first such measurement of a merging cluster was 



2 



Clowe et al. 



performed on 1E0657-56 (aka "The Bullet Cluster"), 
where the gravitational potential of the two merging 
components were found to be spatially coincident with 
the cluster galaxies but significanctly displaced from the 
X-ray plasma clouds (|Clowe. Gonzalez, fc Markevitchl 
12004k iClowe et al.l 12006^ These observations were 
later refined using strong lensing to place tighter con- 
straints on the location and shape of the gravita- 
tiona l potential near th e cluster cores ([Bradac et al.l 
l2006t IBradac et al.1 12001 . and were used to place con- 
straints on the self- interaction cross-section of dark mater 
(jMarkevitch et al.l 12004b iRandall et alll2008f ) as well as 
on any hypothetical "5th force" that affects only dar k 
matter ([Springel fe Farrarl l2007t iFarrar fc Rosen! I2007D . 
The Bullet Cluster has also been used to test alterna- 
tive theories of gravity, with the common result being 
that modified gravity models can reduce the amount of 
dark matter needed in the system, but that the major- 
ity of the mass must still be in a dark, relatively non- 
collis i onal form (lAngus et al.ll2007t iBrownstein fe Moffatl 
120071: iFeix. Fedeli. fc Bartelmannl 12008ft . Several other 
merging clusters have since been found that provide sim- 
ilar results to those from the Bu llet Cluster. Thes e 
clusters are MACS J00254-12 2 2 (IBradac et al~[ 2008) 



A1758 (lOkabe fc Umetsul [2008; Rag ozzine et al 



2012) 



A2744 (IMerten et al.ll201lD . A2163 (ISoucailll2012D. A754 



A175 0, A1914, A2034, and A2142 (jOkabe fc Umetsul 
12001 . 

A weak lensing study (Mahda vi et al.1 120071 

hereafter M07) of A bell 520 (hereafter A520, 

lAbell. Corwin. fc Olowir] [1989), at z = 0.199 
(|Struble fc Roodl [1999) , finds instead a large weak 
lensing signal in the location of the primary X-ray 
plasma cloud (location 3 in Fig. [2]), well away from 
any large concentrations of cluster galaxies, and no 
significant weak lensing signal in the location of one of 
the cluster galaxy concentrations, lab eled as structure 5 
in Fig . [21 A mass reconstruction by lOkabe fc Umetsul 
(2008) using part of the same imaging set as M07 did find 
significant mass in structure 5, but also found sufficient 
mass in location 3 to be consistent with the M07 results. 
More recently, t his "dark peak" has been confirmed by 
Uee et all (|2012l her after J12) using a single passband 
imaging mosaic from the Wide Field Planetary Camera 
2 (WFPC2) on the Hubble Space Telescope (HST). The 
confirmation is not wholly independent as the WFPC2 
mosaic data was combined with a ground based image 
from the Canada France Hawaii Telescope (CFHT) that 
was used in M07. The majority of the signal on the 
locations and masses of the detected structures in the 
core of A520 does, however, come from the new WFPC2 
data in the combined data set. The confirmation of 
the dark peak in J12 is claimed at a ~ 10cr detection 
level, with similar masses for the core structures as was 
measured in M07 except for a detection of significant 
mass at the cluster galaxy concentration of l ocatio n 
5, similar to the results of lOkabe fc Umetsul (2008). 
M0 7 and J 12 provide several scenarios for how such a 
dark peak could arise, such as a filamentary structure 
extending from the cluster, ejection of brig ht galaxies 
from a core during the merger process (e.g. iSales et"aLl 
2007), or a large self-interaction cross-section for dark 
matter, so large that their quoted value lies beyond 
the 5-<7 upper limit on the cross-section derived from 



the Bullet Cluster ([Randall et al.l 120081 ). These results 
have also been used by several aut hors to argue in favor 
of an alternativ e gravity model (|Moffat fc Tothl 120091; 
IBekensteir]|2010D . 

This is not the first time, however, that an appar- 
ently high-significance mass over-density that is not near 
a galaxy over-density has been found in weak lensing 
mass reconstructions . In the cluster A19 42, a mass over- 
density was found bv lErben et all (2000) roughly 7' away 
from the cluster core and not near any significant galaxy 
concentrations. In a separate case, [Miralles et al. (2002) 
found in a blank field STIS image a set of 11 galaxies 
arranged in a pattern reminiscent of those seen in cores 
of massive clusters with multiple strongly lensed back- 
ground galaxies. In both cases, the significance of the de- 
tections, as measured by the likelihood that a randomly 
chosen set of galaxies within the survey area would have 
a similar correlation in their orientations on the sky, were 
sufficiently large that one would not expect to find such 
systems by chance. However, in both cases, deeper ob- 
servations resulted in the measurement of fainter galaxies 
that do not have the same correlated orientation, and 
therefore in the new analyses the mass over-densities 
either greatly diminished in ampl itude (in the case of 
A1942 fvon der Linden et all 120061) ) or co mpletely van- 
ished (in the case of the STIS dark lens ([Erben et al.l 

Eool)). 

We present fully independent weak lensing observa- 
tions of the merging cluster system A520 from a com- 
bined imagining data set from the Magellan 6.5m tele- 
scope in Chile and the Advanced Camera for Surveys 
(ACS) on HST. These ground-based images are of much 
longer exposure times than those used in any of the pre- 
vious weak lensing studies on this cluster and the ACS 
mosaic is deeper than the WFPC2 mosaic due to the 
higher throughput of ACS, and has three observed pass- 
bands for color selection of galaxies as compared to the 
monochromatic WFPC2 mosaic. We investigate whether 
this deeper data set confirms the existence of a significant 
mass over-density at the location of the X-ray plasma 
cloud. The observations are presented in [J2]and the weak 
lensing analysis in $3] Discussion of the results are pre- 
sented in fJH and we summarize our conclusions in ij5j 
Throughout this paper we assume a cosmological model 
with Q m = 0.27, Q A = 0.73, H = 70 km/s/Mpc and 
conventional gravity unless stated otherwise. 

2. OBSERVATIONS 

2.1. HST ACS Images 

We obtained HST imaging with ACS on Feb 25-26 and 
Apr 6-7 2011 (HST Cycle 18 proposal 12253, PI Clowe). 
The new ACS data consist of four pointings in F435W, 
F606W, and F814W. The corresponding exposure times 
are 2300s (1 orbit) in F435W and F606W, and 4600s (2 
orbits) in F814W per pointing. Each orbit was split into 
4 dither positions, with a large enough offset to cover the 
chip gap. 

Since the primary goal of this program is weak lensing 
analysis we took special care when reducing and com- 
bining the images. During its ~ 10 years above the pro- 
tection of the Earth's atmosphere, ACS has accumulated 
significant radiation damage that has degraded its CCD 
detectors. After each exposure, as photoelectrons are 
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transferred through the silicon substrate to the readout 
electronics, a certain fraction is temporarily retained by 
lattice defects created by the radiation d amage, and re- 
leased after a short delay (jJanes ick 2001). This effect is 
known as 'Charge Transfer Inefficiency' (CTI) and spuri- 
ously elongates the shapes of (in particular) faint galaxies 
in a way that mimics weak gravitational lensing. 

We have extended the CTI measurements of IMassev! 
(2010) with a new analysis of hot pixels in extragalactic 
archival HST imaging taken before and after the A520 
data. Interpolating to the two epochs during which A520 
data were obtained, this analysis suggests that ACS ob- 
servations taken in Feb 2011 (Apr 2011), 3283 (3324) 
days after launch suffer from 1.36747 (1.37949) traps per 
pixel, and that the residual effective trap density after 
correction is lower by a factor 20. There was insufficient 
data taken near that time to directly measure the spu- 
rious shear in images obtained at that time. However, 
extrapolating from the trap densities measured in HST 
COSMOS imaging (jMassev et al.ll2010D . this corresponds 
to a spurious shear before correction of ~ 14% spurious 
shear for faint (26 < mp8i4w < 27) galaxies furthest 
from the readout register. This drops rapidly to ~ 4.5% 
by 25 < TOF814W < 26 and further at brighter magni- 
tudes, and falls linearly to zero as one approaches the 
readout register. 

Two independent pipelines have been developed to cor- 
rect the image trailing. Both use the iterative scheme of 
iBristowl J2003) to move electrons back, pix el-by-pixel, to 
where they belong. The first pipeline, by IMassev et al. 
( 20101) (and updated for post-SM4 operations by IMassev 
(2010)) is based around a physical model of charge cap- 
ture and r elease dShocklev fc Readlll952t 15151119521); the 
second, by lAnderson fc Bedinl (|2011l ). is built empirically 
from the observed trail profiles. Both methods have a 
demonstrated level of correction that leaves sub-percent 
spurious shear residuals everywhere on the image at all 
magnitudes. We separately apply each of these pipelines 
to the ACS imaging as the first step in data reduction. 
Furthermore all images taken with the ACS/WFC af- 
ter Servicing Mission 4 show a row-correlated noise due 
to the CCD Electronics Box Replacement . We correct 
for it using the pyraf task acs_destripe (|Grogin et al.l 
12010D FI. The resulting weak lensing shear measure- 
ments for the two CTI correction pipelines were consis- 
tent within 1% of each other, in agreement with our es- 
timate of the expected residual shear from the CTI cor- 
rection. At this level there is minimal effect on the weak 
lensing mass measurements presented herein. 

To stack the correcte d data we use the Multidrizzle 
(|Koekemoer et ail [2002) routine to align and combine 
the images. To register the images we determine the off- 
sets among the individual exposures by extracting high 
S/N objects in the individu al, distortion corrected e xpo- 
sures. We use SExtractor (jBertin fc Arnoutslll996D and 
the IRAF routine geomap to identify the objects and 
calculate the residual shifts and rotations of individual 
exposures, which were then fed back into Multidrizzle. 
We use square as the final drizzling kernel and an output 
pixel scale of 0.05". The resulting images have 5er limit- 
ing magnitudes for galaxies, based on where the number 
counts depart from an exponential growth function, of 

1 http: //www. stsci . edu/hst/acs/sof tware/destripe/ 



TIF435W = 27.0, mp606W = 26.8, and mF8i4w = 26.5. 
2.2. Magellan Optical Images 

We oboserved A520 with the IMACS camera on the 
Magellan Baade telescope during January 16-19, 2004. 
The camera was in the f/4 setup, resulting in a 0V 111 
arcsecond/pixel plate scale and 15.'4 field of view. Dur- 
ing this time there were two significant problems with 
IMACS: The atmospheric distortion corrector had not 
yet been delivered to the telescope, and a problem with 
the CCD amplifiers created horizontal streaking in im- 
ages after a saturated pixel was read. The lack of the 
ADC caused the flat part of the focal plane to be much 
smaller than it was supposed to be, which resulted in 
only the central ~ 6' being in focus, and the image get- 
ting further out of focus the further one moves away from 
the center of the camera. As a result, while many of the 
images were obtained with ~ 0'.'6 seeing, they had l'.'O 
effective seeing at the edges of the image, ~ 8' from the 
center, and l'/4 seeing with noticeable coma in the cor- 
ners of the images. The horizontal streaking occurred in 
4 of the 8 CCD chips, but only after highly saturated pix- 
els were read out. The magnitude of the streaks seemed 
to be dependent on both the total charge in the saturated 
pixels and, oddly, the vertical position of the saturated 
pixel on the CCD - in two of the chips, the streaking 
was very strong at the top and bottom of the chip, but 
almost entire gone in the middle, despite being clearly 
caused by either the horizontal read-out register or the 
on-chip amplifier. Due in part to our limited data set 
and the large changes in the streak amplitude with both 
total charge and chip position, we were unable to find 
a good method of subtracting the horizontal streaking. 
We therefore left it in the images, being sure to mask any 
streaks prior to the creation of flat fields and removed any 
galaxies that overlapped a streak from our weak lensing 
galaxy catalog. 

We observed A520 in three passbands, Bessel B, V, 
and R, with single image exposure times of 5 minutes, 
chosen as a compromise between minimizing the number 
of saturated stars on the images and minimizing the time 
lost to CCD readout. Between each image, we moved the 
telescope by 15" to fill in chip gaps and sample around 
bad pixels. Our final integration times were 120 minutes 
in R and 40 minutes each in B and V. Seeing varied 
between 0"5 and 0"7 in the R images, 0"7 - 0"9 in the 
V images, and ~ l'.'O in the B images. Conditions were 
largely photometric, with stellar fluxes varying by only 
a few percent from image to image. 

We performed image reduction by following the 
prescription for mosaic CCD reduction given in 
iClowe fc~Sc hneider (2001), doing bias subtraction with 
master bias frames, corrections by fitting the overscan 
strip on each chip, and create flat fields by averaging 
together the science images with sigma-clipping after re- 
moving all detected objects from the images. We register 
the images using a two step process of converting each 
CCD to a detector plane coordinate grid using a linear 
shift in x and y and a rotation in the x-y plane. We 
then map the detector plane onto the sky using a 7th 
order two-dimensional polynomial by com paring stellar 
posit ions to those in the USNO-B catalog (|Monet et al.l 
2003). All of the images use the same CCD to detec- 
tor plane conversion parameters, but the coefficients of 
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the detector plane to sky conversion polynomial freely 
vary for each image. We therefore have 21 free param- 
eters from the CCD to detector plane conversion and 
36 x JT-images free parameters from the polynomial coef- 
ficients, and roughly 200 x ni magcs stellar positions to 
constrain the fit. The resulting stellar positions have an 
average root mean squared (rms) of 0'.'004 compared to 
the same stars in other images from this dataset, and 
a rms position difference of 0'.'25 when compared to the 
USNO, which is fairly typical of t he positional uncer - 
tainties within the USNO-B catalog ([Monet et al.l f2003). 
One possible source of failure in this method is if the 
CCD chips in the camera are not sufficiently well aligned 
vertically, because then the detector-plane to sky coor- 
dinate conversion can change too rapidly across the chip 
gap for the relatively low-order polynomial. We test for 
this by comparing the rms positions of stars that appear 
on more than one chip to those residing exclusively on 
a single chip and by looking for changes in the shapes 
of the PSF across chip gaps. In both tests, we find no 
significant deviation that would indicate a vertical mis- 
alignment of the CCDs to a degree that would affect 
either the image registration process or the subsequent 
weak lensing analysis. 

Using the polynomials from the registration process, 
we map the images onto a common coordinate grid, pre- 
serving the 0'/ 111 arcsecond/pixel plate scale and ori- 
entation of the original images by using a triangular 
method with linear interpolation that preserves surface 
brightness and has been shown to not induce system- 
atic changes in ob ject shapes for fractional pixel shifts 
([Clowe et al.l 12000ft . We produce the final images by co- 
adding the registered images using a sigma-clipping algo- 
rithm to detect and remove cosmic rays, while not clip- 
ping the centers or wings of stars. The final images have 
FWHM in the central 6' of Of! 63 in R, 0'/75 in V, and 
l'/05 in B, with increasing FWHM with distance from 
the center of the image. The 5<r limiting magnitudes 
for galaxies in the image centers in the final images, as 
measured from where the number counts depart from an 
exponential growth function, are mj = 26.2, my = 25.8, 
and iur = 25.7. The 5<r limiting magnitudes for galaxies 
at the edges of the images decrease by ~ 0.2 magnitudes 
due to the larger PSF size. 

2.3. X-ray Images 

We have created an approximate projected gas mass 
map using a 0.8-4 keV X-ray image extracted from the 
archival Chandra 520 ks dataset (Markevitch et al. in 
preparation). The X-ray emissivity at photon energies 
E <C T e depends very weakly on gas temperature and 
its variations across the cluster. A520 is a T e — 7 keV 
cluster, and Chandra has a peak of sensitivity at E ~ 1 
keV, which makes the X-ray surface brightness in our 
energy band a good representation of the projected X- 
ray emission measure, EM oc n e n p . 

To convert the projected emission measure to the gas 
mass requires knowledge of the three-d imensional cluster 
geometry. Unlike the B ullet cluster ([Markevitch et all 
l2002t iClowe et al.l [2006) that appears to have a simple 
geometry, A520 is irregular and we cannot make any 
plausible assumptions about its gas distribution along 
the line of sight. An approach often used in such sit- 
uations to obtain a first-approximation gas mass map 



is to take a squar e root of the X-ray brightness. As 
in Ragozz ine et al.l (|2012[ ). we attempt a slightly higher 
level of accuracy by taking advantage of the fact that 
clusters are centrally peaked and approximately spher- 
ically symmetric on large scales. To do this, we first 
fit a spherically-symmetric /3-model to the X-ray radial 
brightness profile and create a projected gas mass that 
corresponds to that model. This zero-approximation 
mass map is then multiplied by a factor (Sx/Sp) 1 / 2 , 
where Sx is the cluster surface brightness and Sp is the 
/3-model image. To nor malize this gas mass, w e use a 
A/gas - T relation from IVikhlinin etall (|2009ft . which 
was derived from the Chandra gas masses and X-ray 
temperatures, and the overall cluster temperature of 7.1 
keV (Govoni et al. 2004). Though the M gas - T rela- 
tion is derived for relaxed clusters, hydrodynamic simu- 
lations indica te that it should not be very differe nt for 
mergers (e.g. INaeai. Vikhlinin. fe Kravtsovl 12007). For 
the A520 temperature and redshift, the relation gives 
-Mgas = 7.9 x 10 13 M© in a sphere of radius r 5Q0 = 1.10 
Mpc. We normalize our map to have the same gas mass 
within the r$oo aperture as that for the /3-model. For 
a check, we have also tried a more direct (but also more 
noisy for such irregular clusters as A520) estimate for 
the normalization for our /3-model using the A520 Chan- 
dra spectrum from the central r = 3' region and fittin g 
it with the APEC spectral model (|Foster et all 12012ft . 
whose normalization gives the projected X-ray emission 
measure. This gives a gas mass 18% higher than the 
above value within the same sphere. 

Simulations, e.g., by Kravtsov et al. (2006) and Rasia 
et al. (2011), indicate that even the extreme merging 
clusters, such as A520, follow the M — T relation with 
a scatter of about 20-25% along the mass axis, and the 
Mgas — Mtot relation is even tighter. Other errors in 
our analysis should be smaller, and we have assigned 
a conservative 25% error (68% confidence) to the gas 
masses. The above 18% difference between our two gas 
mass estimates is well within this assumed uncertainty. 

3. WEAK LENSING ANALYSIS 

3.1. Shear Measurement 

We perform weak lensing analysis on the images with 
the goal of obtaining a two-dimensional distribution of 
the surface density in the cluster. This is done by 
measuring second moments of the surface brightness to 
calculate an ellipticity for each galaxy, correcting this 
shape for smearing by the point spread function (PSF) 
to measure a shear, rejecting stars by size and central 
surface brightness, and rejecting likely cluster and fore- 
ground galaxies by color. The methodology we use for 
the PSF correction is that of a modi fied KSB technique 
([Kaiser. Squires, fc Broadhurstl 119951) . details for which 
can be found in IClowe et al. 1 120061) . Galaxies selected 
for the weak lensing analysis had a photometric S/N > 10 
in the R or F814W passbands, R > 22 or F8UW > 21.5, 
and did not have any bright neighbors near enough to 
significantly influence the second moment measurements 
(less than 3 times the sum of the scale radii of the galaxy 
and the neighbor). Weak lensing analysis is performed 
separately on the co-added Magellan image and on each 
of the four ACS pointings, then we combine the result- 
ing shear catalogs to produce a final catalog. The weak 
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Figure 1. Shown above in the left (right) panel is the distribution in color-color space of galaxies in the Magellan (HST) images. The 
cyan points show all of the galaxies with magnitude R < 24 (F606W < 24), while the black, magenta, and blue asterisks denote the colors 
of spectroscopically confirmed cluster, foreground, and background galaxies respectively. The red circles show colors of stars, rejected from 
the galaxy catalogs using size and central surface brightness criteria. The solid box shows the selection for likely cluster galaxies used to 
perform the cluster luminosity measurements, while the dashed lines show the cuts used to remove likely cluster and foreground galaxies 
from the shear catalog (galaxies to the right and below the dashed line are kept). The red and bl ue lines plotted on the right panel show 
theoretical colors for galaxies formed at z = 6 with Salpeter initial mass function (Salpctcr 1955) and present day solar metallicity with 
a single star-burst population (red) and an 10 Gyr exponential decay star formation rate (blue). The tick marks and associated labels 
indicate where galaxies of a given redshift reside along these lines. The color evolution models were generated using the EzGal software 
(Manconc & Gonzalez 2012) using the CB07 models. 



lensing measurements were performed with a modified 
version of the IMCAlQ software package. 

For the Magellan image, two additional defects in the 
image require added modification of the shear measure- 
ments. The first is the horizontal streaking coming off of 
saturated stars in 4 of the CCDs. As these streaks ap- 
pear to be a change in the bias and/or gain of the readout 
amplifier, we do not trust the shapes of any objects in 
these regions, and therefore simply remove all objects in- 
tersecting any of these horizontal streaks in the image. 
Due to the non-local nature of gravitational shear, this 
removal of galaxies will not bias our results, except to 
slightly increase the noise in the mass reconstructions in 
the vicinity of the removed galaxies. 

The second defect is that the images go out of focus 
at the edges of the image, which causes a strong change 
in the shape and size of the PSF as a function of radial 
position from the center of the image. For the i?-band 
image, from which we measure the galaxy shapes, the 
PSF size (FWHM) increases from ~ 0'.'63 in the center 
to ~ 0'.'85 at the edges and ~ l'/4 in the corners. Fur- 
ther, coma can be seen in the PSF in the corners of the 
image. As a result, we restrict our shear measurements 
to a 8' radius from the center of the image and, instead 
of a straight size cut to separate stars and galaxies, we 
use a 7th order polynomial fit to the stellar half-light 
radius as a function of image position and rejected any 
object with a size smaller than O'.'l larger than the fit 
value at that location. We supplement the stellar rejec- 
tion by also rejecting objects with unusually high central 
surface brightness for its magnitude. We also measure 
the KSB PSF correction terms (P s h,-Psm, and stellar el- 
lipticity) using a broad range of weighting function sizes, 

2 http: //www. if a.hawaii . edu/kaiser/imcat 



fit these as 7th order polynomials for image position vari- 
ations, and use the fitted values for the PSF correction of 
a given galaxy based on its position and size. To obtain 
the final P 1 correction factor in the KSB technique, we 
divide the image up into four regions based on PSF size 
to fit P 1 as a function of galaxy size and ellipticity and 
reduce the significant noise present in the P 7 measure- 
ment for each individual galaxy. Using simulations with 
PSFs taken from the Magellan image, we found this tech- 
nique systematically underestimates the measured shears 
by - 13% for a 0'.'6 PSF increasing to - 15% for a 0"8 
PSF for the smallest galaxies in the simulations, decreas- 
ing to ~ 10% for galaxies significantly larger than the 
PSF size. We determine this correction factor from fits 
to the simulation results based on the size of the PSF in 
the galaxy's location and the size of the galaxy. 

For the HST ACS images, we rejected stars using a 
size cut (< 0'.'0.081 for 50% encircled light radius) as well 
as rejecting objects with unusually high central surface 
brightness for their magnitude. We again measure the 
KSB PSF correction terms for a range of weighting 
function sizes, fitted these using a 5th order polynomial 
for image position variations in each pointing, and use 
the fitted values matched to the galaxy size for correct- 
ing the PSF smearing. We use the ACS- like STEP3 
( |http://www.roe.ac.uk/^heymans/step/cosmic_shear_test.htm.i| 
simulations to calibrate the PSP corrections, finding 
a systematic underestimate of ~ 8% for the shear 
measurements, which was corrected for in the ACS 
measurements. We perform the shear measurements 
independently for each of the three ACS passbands. 

Because the HST and Magellan images observe galaxy 
populations with different redshift distributions, due 
mainly to loss of shape information on intrinsically 
smaller galaxies with the ground based PSF, and the 
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strength of the shear measurements depends both on the 
mass of the l ens and the redshifts of the background 
galaxies (see ^3.2[) . before the shear measurements can 
be averaged between the two datasets we need to adjust 
the catalogs to have the same mean lensing depth. Using 
external photometric redshift catalogs (see §3.41 for de- 
tails), we determined that the HST dataset would have 
a mean lensing signal in a given region that is ~ 1.05 
times that of the Magellan image. We therefore scale the 
shears measured in the Magellan image by 1.05 before 
combining with the HST catalog to create the final weak 
lensing catalog. 

We compute weights for each galaxy in each data set 
by computing the inverse of the rms shear for nearby 
neighbors in significance and size space, with each data 
set showing that large, bright galaxies have a rms in- 
trinsic shape of rms s =0.24 per shear component in the 
F814W ACS passband, 0.26 for the F606W ACS pass- 
band and the Magellan i?-band, and 0.27 for the F435W 
ACS passband. Fainter and smaller galaxies have larger 
rms shear values, indicating increasing measurement er- 
rors for the second moments from sky noise and PSF 
correction factors. We therefore separate the rms shear 
values into two components, an intrinsic shape value cho- 
sen to be rmSi n = 0.24 and a measurement value com- 
puted as rms m = ^/rms^ — rms? n , and set a lower limit 

on rms m = 0.05 based on the spread in the rms 9 val- 
ues for the brightest and largest galaxies. From these we 
create two weighting functions, w g — l/rms g for weight- 
ing shears in the weak lensing mass reconstructions, and 
w rn = l/rms m for weighting the co-addition of shears 
for galaxies with multiple shear values in different ACS 
passbands, overlapping ACS pointings, and those galax- 
ies located in both the ACS and Magellan images. We 
compute a final weight for each galaxy by adding the w m 
values for each shear measurement in quadrature, taking 
the inverse to get a final rms m , and taking the inverse of 
rms m added in quadrature with rmsi n . 

The final step in creating the weak lensing shear cat- 
alogs is to remove likely foreground and cluster galaxies 
from the galaxy catalogs. T his we do by using the Hyperz 
photo metric redshift code (jBolzonella. Miralles. fc Pellol 
2000) to produce magnitudes in each observed passband 
for a range of galaxy templates from starbursts to passive 
ellipticals for < z < 0.25, adjusting these for the galac- 
tic extinction of the A520 field, isolating the regions in 
color-color space for the ACS and Magellan passbands 
for these galaxies, and removing all galaxies from the 
shear catalog that have colors, within photometric er- 
rors, which lie within the low-redshift galaxy color-color 
regions. To account for noise in the photometric mea- 
surements of fainter galaxies possibly moving foreground 
and cluster galaxies across the selection boundaries, we 
excluded all galaxies whose colors were within la of the 
boundaries (a more stringent cut at higher a resulted in 
the loss of too many faint galaxies and a severe decrease 
in the S/N of the shear measurement). These color-color 
cuts, shown in Fig. [I] remove ~ 40% of the Magellan and 
~ 30% of the ACS galaxies that otherwise were consid- 
ered bright enough, large enough, and isolated enough 
to provide good shear measurements. The final lensing 
catalog has a number density of galaxies of 22 per square 
arcsecond for the Magellan images and 56 for the ACS 



images, which result in rms shear per square arcminute 
of 0.036 and 0.058 for the regions around the core of 
A520 with and without ACS imaging respectively. The 
number density of galaxies decreases by ~ 15%, and rms 
shear increases by ~ 10%, in the Magellan image as one 
approaches the edges of the image due to the increased 
PSF size. 

3.2. Mass Reconstruction 

The PSF corrected galaxy ellipticity measurements 
each provide independent measurements of the re- 
duced shear q, where q = 7/(1 — «) (see 
iBartelmann fc Schneider] 120011 for formal derivations 
of the weak lensing concepts). The shear 7 is the 
anisotropic change in the light distribution of the back- 
ground galaxy caused by weak lensing, and the conver- 
gence k is the isotropic change in the background galaxy's 
light distribution. The convergence is also the surface 
density of the lens (£) scaled by a geometric factor (S C rit) 
that depends on the angular diameter distances between 
the observer, the lens, and the background galaxy being 
lensed. Thus, to study the mass distribution of A520 
with weak lensing, we need to convert our measured re- 
duced shear data points to the convergence field of the 
lens, use an assumed cosmology and mean redshift for 
the background galaxies to estimate £ C rit, and scale the 
convergence field to a measurement of the surface density 
of the lens. 

To convert the reduced shear measurements to the 
convergence field, we u sed the iterative technique of 
iSeitz fc Schneider] (119951). This technique is based on 
the KS93 algorithm (|Kaiser fc Sauireslll993[ ). which uses 
that both 7 and K are second derivatives of the surface 
potential to combine derivatives of 7 to get derivatives of 
K, which are then integrated to produce the convergence 
field to within an unknown integration constant, which is 
the mean value of k at the edge of the reconstructed field. 
Because the input galaxy catalogs provide only a sparse 
sampling of the reduced shear field, the output k field 
needs to be smoothed to remove large noise spikes, in 
this case by convolution with a a = 60kpc Gaussian ker- 
nel. The iterative technique is to initially assume k = 
across the reconstucted area, so 7 = g, obtain a measure- 
ments of the k field, and then use this to perform a new 
correction of 7 = g X (1 — k). After four iterations, we 
find the difference between the input n field and the out- 
put K field from the KS93 algorithm differ by less than 
0.01% of the input field, and stop the iteration. For our 
combined catalog of 5903 background galaxies, the full 
set of iterations to produce a final two-dimensional (2-D) 
mass reconstruction takes only a few seconds. 

The integration constant in each reconstruction is cho- 
sen by letting the mean value of k at the edge of the 
reconstructed area equal to the expected density of a 
cluster with the observed X-ray temperature. For our 
A520 data, this sets the convergence of the lens at a ra- 
dius r ~ 1500 — 1600 kpc to n ~ 0.01, which is typical 
of clusters with X-ray te mperatures of ~ 7 keV b ased on 
the Af 500 - T x relation of iVikhlinin et al.l (|2009f >. While 
there is some dependence on the mass measured in the 
core region of A520 on this outer value, varying the con- 
vergence at this outer radius by ±0.01 results in a change 
in our measured substructure masses of less than 5% ab- 
solute mass and 2% relative mass ratios of structures, 
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Figure 2. Shown above in the left panel is a 7' X 6.'75 color composite from the HST ACS mosaic images with the weak lcnsing surface 
density reconstruction overlayed in cyan contours and the Chandra X-ray derived gas surface density in magenta contours. The weak 
lensing contours show steps in surface density of 2 X 10 s M0/kpc 2 (k steps of 0.056) above the mean surface mass density at the edge 
of the Magellan image (~ 1600 kpc radius), and the gas mass contours show steps of 7.4 X 10 6 M©/kpc 2 with the outer contour starting 
at 4.4 X 10 7 M Q /kpc 2 . The upper right panel shows the weak lensing contours superimposed on a smoothed cluster galaxy luminosity 
distribution in greyscale, with both the luminosity and surface density distributions smoothed by the same a = 60 kpc Gaussian kernel. The 
bottom right panel shows contours of the mass aperture statistic from the weak lensing data, with contours of steps of lcr, superimposed 
on the cluster galaxy luminosity distribution. Also labeled in the left panel are the regions of structures 1-6 identified in M07 and J12 as 
well as the new structure 7. 



and has no discernible impact on the shapes or centroids 
of the substructure mass peaks. 

The resulting k distribution is shown in contours over- 
layed on a color image constructed from the ACS mosaic 
data in the left-hand side of Fig. [2] In the upper right- 
hand panel, the k contours are overlayed on a greyscale 
map of the luminosity distribution of cluster galaxies, se- 
lected by using the same color-color cuts that were used 
to exclude likely cluster galaxies from the weak lensing 
galaxy catalog, smoothed by the same sized Gaussian 
kernal as is the mass reconstruction. We detect four pri- 
mary mass concentrations (1, 2, 4, and 5 in Fig. [2]), and 
see some evidence for excess mass in region 6 although 
none of our later tests would argue for a significant detec- 
tion of an additional cluster substructure in this position. 
We find no evidence of the mass overdensity in the dark 
peak region 3, and instead find a surface density distri- 
bution in that region which is in very good agreement 
with the underlying cluster galaxy light distribution. 

In the bottom-right panel of Fig. [21 we also s how the re- 
sults of using a mass aperture (M ap ) statistic (jSchneiderl 
1996) on the weak lensing catalog. The M ap statistic 
measures the k distribution convolved with a compen- 
sated filter, and has an advantage over the 2-D mass re- 
construction in that it produces easily measurable errors. 
The disadvantage is that the M ap statistic uses a more 
limited radial extent of the reduced shear measurements, 
and therefore has a lower signal-to-noise (hereafter S /N) 
in its measurement than the 2-D mass reconstruction. 



To avoid having the negative portion of the compensated 
filter overlapping nearby structures, and thereby signifi- 
cantly decreasing the M ap signal, we used a 200 kpc outer 
radius for the statistic, and measured the values for cen- 
ters distributed on a 100 x 100 grid across the cluster 
core region shown in Fig. [2j We detect at > 2a signifi- 
cance structures 1, 2, 4, and 5, detect structure 6 at only 
la due in part due to its proximity to structure 4, and 
again find no detection of excess mass in the vicinity of 
the proposed dark peak. We do find an additional struc- 
ture that we label as structure 7 in Fig. [2j however it 
is significant only in the M ap measurements and not in 
the full mass reconstruction, which would be consistent 
with it being a local noise peak caused by a small num- 
ber of highly elliptical galaxies in the vicinity. We do not 
otherwise consider structure 7 in this paper. 

3.3. Bootstrap Resampling 

Determining the errors on the 2-D mass reconstruction 
is more problematic than the errors on M ap , as varia- 
tions in the number density and magnitude of the intrin- 
sic cllipticity of background galaxies cause the errors in 
both the enclosed mass, the mass centroid, and the mass 
structure shape to all vary by large amounts across the 
reconstructed area. A common, but incorrect, method 
used to estimate these errors is to measure the rms shear 
and the mean density of the background galaxy catalog, 
and propagate these errors through the mass reconstruc- 
tion algorithm obtaining an average noise level for the 
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Figure 3. Shown above are nine randomly-selected weak lensing surface density reconstructions from the 100,000 bootstrap resampled 
catalogs used to measure errors in the weak lensing reconstructions superimposed on the cluster galaxy luminosity distribution. The 
contour levels are identical to those in Fig. [2] and the magenta x shows the location of the dark peak of J12 (structure 3). The middle right 
reconstruction shows a structure that is morphologically similar to and has similar mass as the dark peak in J12; such structures are found 
in ~ 2% of the bootstrap resampled reconstructions, and are the only reconstructions in which the column mass in the dark peak location 
agrees with that of J12. Other reconstructions show smaller peaks near the dark peak, but have much less mass than reported by J12. 



reconstruction. The problem with this approach is that 
the k measurement for a given peak location is measured 
from the shear of galaxies in the catalog with an effec- 
tive weighting of j/r. As 7 is largest and r is smallest 
for galaxies immediately around the peak location, most 
of the weight in the n determination comes from a rel- 
atively small number of the nearest galaxies. If any of 
these galaxies have an intrinsic ellipticity near the edges 
of the distribution function, the noise in that peak will 
be significantly larger than average. 

Another method to generate random noise fields that 
is commonly used is to preserve the position and total 
magnitude of the reduced shear measurement for each 
galaxy, but to apply a random orientation to each galaxy 
before performing the mass reconstruction. After doing 
this enough times, one can then compute a rms of the 
noise field in each pixel of the reconstruction. The prob- 
lem with doing this in a field containing a massive cluster 
is two-fold: the cluster shear is still part of the measured 
reduced shear, so one would be significantly overestimat- 



ing the combined intrinsic shape and measurement noise 
for galaxies near the cluster core, and the average k in 
the random reconstructions will be 0, and thus the re- 
construction will be misinterpreting the reduced shear 
to be shear, which also artificially enhances the level 
of the noise in the vicinity of the cluster core. Thus, 
this method will give a lower limit on the significance 
of the detection of a structure in the weak lensing mass 
reconstruction, but in simulations of massive clusters we 
often find the lower limit can be as low as half of the 
true significance. One can try to correct for this by us- 
ing the smoothed mass reconstruction from the data to 
change the reduced shear into shear and then subtract 
off a shear field created from the mass reconstruction. 
By doing this, however, you artificially reduce the level 
of the noise in the vicinity of noise peaks in the mass re- 
construction, and thus overestimate the significance, in 
simulations often as much as a factor of 2, of the noise 
peaks. One can therefore use this method to measure a 
minimum and maximum significance for the structures 
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Figure 4. Shown above are contours enclosing the mass centroid 
position uncertainties for structures 1 (red), 2 (blue), 3(magenta), 
4 (green), and 5 (cyan) superimposed on the cluster galaxy lu- 
minosity in greyscale. The contours enclose the locations of the 
detected mass centroids in 100,000 bootstrap resamplings of the 
reduced shear catalog, and encompass 68%, 95%, and 99.7% of the 
centroid measurements. 

in the reconstruction, but the range between these two 
estimates is often quite large. 

To properly measure how noise in the weak lensing 
catalogs affects the 2-D reconstruction, one needs to use 
a method that preserves the underlying reduced shear 
field while simulating the noise, which comes primarily 
from the intrinsic ellipticity of the background galaxies. 
The method that we use to do this calculation is boot- 
strap resampling of the background galaxy catalog, in 
which one creates a new catalog with the same number 
of entries as in the original catalog, with each entry be- 
ing a randomly selected member of the original catalogs 
and objects are allowed to be selected more than once. 
For a suitably large catalog, this results in the chance of 
any given object having an integer weight m > to be 
e~ 1 /ml, and the chances of any group of n galaxies not 
being in the new catalog to be e~ n . Once the new cata- 
log has been generated, a 2-D mass reconstruction can be 
measured from it, and by repeating the resampling as of- 
ten as necessary, suitable statistics on the enclosed mass, 
mass centroid, and structure shapes can be measured. 

We show nine randomly selected bootstrap resampled 
mass reconstructions in Fig. [3] as contours superimposed 
on the cluster luminosity distribution in greyscale. As 
would be expected, in general, the larger the mass peak 
in the original reconstruction, the smaller the changes 
in the relaltive size, position, and shape in the boost- 
rap resampling reconstructions. We see that most of the 
movements in the peak locations are well correlated with 
the shapes of the underlying galaxy distributions, consis- 
tent with a model where the galaxies are tracers of the 
dark matter mass distribution, and the observed peak lo- 
cations are simply the largest noise peak in the vicinity 
of the structure core. In most of the resamplings, we do 
not see a dark peak in the vicinity of location 3, although 
in about 2% of the cases we do see a structure which one 
would identify as a mass peak not associated with cluster 
galaxies between locations 2 and 4. One such example 



can be seen in the middle-right reconstruction of Fig. [31 
Because the mass reconstructions have been smoothed 
to eliminate noise spikes from the reduced shear field, 
using a peak finder to detect the locations of the high- 
est values of k within the various structures is equiva- 
lent to finding the mass centroid of the structure with a 
weighting function equal to the smoothing function used 
in the reconstruction. In Fig. [4] we show the locations 
of the primary mass peaks for the various structures 
in 100,000 bootstrap resamplings, as contours that en- 
close 68%, 95%, and 99.7% of the peak locations, which 
were detected as being the nearest significant peak to 
the structure in the original reconstruction. Structure 

2 is detected in all of the resamplings, and has a major 
axis in its centroid uncertainty distribution that agrees 
well with the NW-SE major axis of the distribution of 
cluster galaxies. Structure 4 is also found in all of the re- 
samplings, with a major axis in its centroid uncertainty 
running mostly N-S, and has a good agreement with the 
location of the cluster galaxies in about 1 /3rd of the re- 
constructions. This result suggests that the northern off- 
set of structure 4 from the galaxy peak seen in the orig- 
inal reconstruction is likely just due to noise and not a 
significant feature of the merger. 

In about 2% of the reconstructions, structure 1 is not 
detected and instead the nearest significant peak is that 
of structure 2, otherwise the location of peak 1 is very 
centralized around the single giant elliptical galaxy in the 
region. Structure 5 shows a major axis in its positional 
uncertainty that runs nearly E-W, in good agreement 
with the distribution of the cluster galaxies in this region. 
The centroid of structure 5 overlays the brighter elliptical 
galaxies on the eastern end of the structure only ~ 10% of 
the time, while being detected in the vicinity of the dark 
peak location 3 about ~ 2% of the time, and not detected 
at all about ~ 1% of the time (when the nearest peak 
is that of structure 2). In the ~ 2% of the cases where 
a peak is found in the vicinity of location 3, it is almost 
always (~ 90%) the case in which peak 5 is located at the 
extreme western end of its positional distribution rather 
than finding separate peaks in locations 3 and 5 as was 
seen in J12. 

We also show in Fig. @] the nearest peak to the dark 
peak location 3 from J12. As opposed to the other 4 
primary structures in the cluster, this distribution has 
multiple peaks in the centroid location distribution and 
is consistent with our finding in each reconstruction the 
nearest noise peak superimposed on a bridge structure 
traced out by the cluster galaxies. When a peak with 
a mass equivalent to that in J12 is found near region 

3 (~ 1% of the reconstructions), 85% of the time it is 
the case that structure 5 is at the extreme western end 
of its distribution and is the primary contributor to this 
increase in mass in region 3. In the bootstrap resampling 
reconstructions, we do not see any indication of a dark 
peak in region 3 that is being suppressed in our original 
reconstruction by the chance projection of a handful of 
highly elliptical objects. 

To determine how these results are influenced by the 
size of the smoothing/weighting function used, we re- 
peated them using 30 kpc and 120 kpc radius smooth- 
ing functions. As expected for a white noise field super- 
imposed on an underlying signal, we detect more mass 
peaks with the 30 kpc smoothing radius than for the 60 
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kpc smoothing, with the location of the most significant 
peak having a larger variation in position than with the 
60 kpc smoothing radius. For the 120 kpc smoothing 
radius, we get a slightly smaller (~ 80 — 90% of the con- 
tour sizes seen in Fig. |4j spread in the centroid locations 
when the structure is detected. Structures 1 and 5 are 
not found ~ 20% of the time, instead appearing as ex- 
tensions of structure 2, and structure 2 is not detected 
~ 5% of the time. Thus, the contours in Fig. [4] are likely 
to be slight overestimates of the true uncertainty in the 
mass centroid locations, but regardless of the smoothing 
functions we never detect a significant peak in the mass 
reconstruction at the location of structure 3 in more than 
5% of the bootstrap reconstructions. 

3.4. Mass Measurements 

Instead of asking whether A520 has a dark peak, a 
more direct question is whether there is excess mass in 
the vicinity of the reported location of the dark peak. 
To do this, we measure a column mass enclosed within 
a given radius by simply integrating over the k values 
within the given radius around a chosen center in the 2-D 
mass reconstructions, and do the same in the bootstrap 
resampled reconstructions to look at the distribution of 
errors in the measurement. For easy comparison, we use 
the same aperture size used in J12, 150 kpc, which results 
in non-overlapping mass integration regions around the 
peaks except for a 9 kpc overlap between peaks 3 and 5 
and a 50 kpc overlap between peaks 4 and 6. The mass 
of the X-ray plasma and integrated cluster luminosity in 
each region are computed using the cutout regions from 
the X-ray mass and cluster light images described earlier. 

To convert the integrated n values to column masses, 
we need to assume a value for S cr ;t. Because we do not 
have data in enough passbands to measure reliable pho- 
tometric redshifts for the background galaxies, we use 
photometric redshifts from other imaging data sets for 
the magnitude range of the background galaxies. From 
these, we calculate the inverse of the mean value of S~| t , 
using weights based on the F814W magnitude to mimic 
the weights used in the weak lensing measurements, that 
such a redshift distribution would have if it were to be 
lensed by a z — 0.2 cluster. The two photometric red- 
shift datasets we u sed were the COSMOS field catalog of 
lllbert et al.l ((2009), for which the photometry is mainly 
ground-based Subaru data and thus a good match to 
the M agellan image, and the UDF catalog of ICoe et al.l 
(2006), which has entirely space based photometry, and 
is thus a good match to the HST mosaic. From these 
catalogs, we determined mean values of £ C rit of 3.4 x 10 9 
M /kpc 2 for the HST mosaic and 3.6 x 10 9 M Q /kpc 2 for 
the Magellan image. As we scaled the Magellan shears 
to compensate for this prior to coaddition of the cata- 
logs, we adopt the UDF value for converting k to surface 
density. The mass ratios of the various apertures are in- 
sensitive to the adopted value of £ C rit, and show only 
minor variations when we change the scaling factor be- 
tween the Magellan and HST catalogs. 

The resulting mass and luminosity measurements for 
the six structures detected in our mass reconstruction 
and the dark peak location can be found in Table 1. 
We also provide the errors for the lensing masses, de- 
termined from the variation in masses measured in the 
bootstrap resampled reconstructions. We do not provide 



errors for the gas mass or cluster luminosity measure- 
ments, as the random errors between the measurements 
for the different locations are dwarfed by the error in the 
weak lensing masses. We also list in Table 1 the mass-to- 
light ratio, computed after subtraction of the X-ray gas 
mass, and its error for each of the peaks. There arc likely 
systematic errors in all three measurements, however the 
systematic errors in the cluster luminosity measurements 
(e.g. galactic dust, intracluster light) will only change the 
absolute scale of the mass-to- light measurements and not 
the relative values of the peaks. Some non-cluster galax- 
ies are likely to have been included by the color selec- 
tion process, however in all 7 regions the majority of the 
galaxy light is emitted by red-sequence galaxies. Using 
the same color selections on the COSMOS and UDF cat- 
alogs discussed above suggests that interlopers should be 
contributing ~3x 10 9 L (calculated by converting flux 
to luminosity by placing all galaxies are at the cluster 
redshift) on average in each aperture, or less than 2% 
of the measured luminosity. Therefore fluctuations of a 
factor of a few in the interloper population among the 
apertures will produce minimal contamination of mass- 
to-light ratios. The systematic errors in the weak lensing 
(e.g. chosen value of S crit , PSF correction) are relatively 
small (< 10%) and while the systematic errors in the gas 
mass measurement could be quite large (~ 25%) the gas 
mass is sufficiently small compared to the weak lensing 
mass in all of the peaks that combined they will cause 
the relative mass-to-light ratios among the peaks to vary 
by less than the random errors. 

From the detected mass structure (1, 2, 4, and 5), we 
obtain a weighted mean mass-to-light ratio of 108 ± 24 
Mq/Lq, which is in excellent agreement with the mass- 
to-light ratio measured for the region of the dark peak 
(114 ± 54 Mq/Lq). For structures 1 — 5 (we leave out 6 
due to the large overlap with the structure 4 aperture), 
the hypothesis that all of the structures share the same 
mass-to-light ratio gives a reduced x 2 = 0.75. If struc- 
ture 4 is assumed to be a noise peak superimposed on 
the underlying mass distribution and instead we measure 
the mass-to-light centered on the cluster galaxy luminos- 
ity peak, the weak lensing mass decreases by only ~ 2% 
while the cluster luminosity increases by ~ 17%, decreas- 
ing its mass-to-light ratio from 133 ± 20 to 111 ± 17, and 
the reduced x 2 °f the constant mass-to-light hypothesis 
to 0.43. None of the other structures are misaligned with 
the underlying light distribution enough to significantly 
change the reduced \ 2 by centering them on their cluster 
luminosity centroids. The mass-to-light ratio of structure 
6 is much higher than the others, but is still within 2a 
of the mean. 

Another method of measuring mass within a given ra- 
dius, and that used by J12, is apertu re densitometry 
(jFahlman et al.lfl99i IClowe et al.l l2000h . The two tra- 
ditional statistics measure the mean n within radius r\ 
minus the mean k within a given annular region, with 
the difference being how the subtractive annular region 
is defined: 



C(ri) = R(r < n) - R(n < r < r max ) 

= -■ _ 2/ 2 — / {lT)d\x\r 
' 1/ 'max J r\ 
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Figure 5. Shown above are the weak lensing surface density reconstructions using only the galaxies detected the Magellan image (left) 
and HST mosaic (right) as contours superimposed on the cluster galaxy luminosity in greyscale. The Magellan reconstruction has been 
smoothed by a 82 kpc Gaussian kernel, while the HST reconstruction was smoothed by a 26 kpc Gaussian kernel. The magenta x shows 
the location of the dark peak in J12. 



Table 1 

Mass Reconstruction Substructure Properties (r < 150 kpc) 



Substructure 


RA 
(h : m : s) 


Dec 
(":': ") 


Column Mass 
(h^lO^Af©) 


Luminosity 
(hTO a 10 n L, @ ) 


Gas Mass 
(hro /2 W 13 M @ ) 


M/L 

(M o /L z0 ) 


PI 


04:54:19.60 


+02:57:49.09 


3.03+0.69 


2.43 


0.25 


114+28 


P2 


04:54:14.84 


+02:57:06.25 


4.08+0.73 


4.16 


0.40 


88+18 


P3 


04:54:11.25 


+02:55:37.28 


2.26+0.75 


1.38 


0.69 


114+54 


P4 


04:54:04.57 


+02:53:58.60 


4.64+0.63 


3.11 


0.50 


133+20 


P5 


04:54:17.11 


+02:55:30.09 


3.00+0.77 


2.66 


0.44 


96+29 


P6 


04:54:09.61 


+02:53:55.90 


3.03+0.66 


1.15 


0.65 


207+57 



Note. — M/L is calculated after subtraction of the gas mass. 



Cc(n) = n(r < ri) - «(r 2 < r < r max ) = 2 1 (7 T )dlnr 

Jn 

+ , _ 2/2 / (lT)d\nr, 

with r max being the maximum radius used in the mea- 
surement, ri is the inner radius of the subtractive aper- 
ture for Cc, 7t is the portion of the shear measurement 
oriented tangential to the chosen center, and the angu- 
lar brackets indicate azimuthal averaging. In practice, 
one can bin the tangential shear measurements in bins of 
constant logarithmic radius change (dlnr) and convert 
the integrals in the above equations to summations over 
the bins. This allows for an easy calculation of the er- 
rors, with a\ = (Ebin 4 ( dln? ') 2 rms2/n bin )/(l-r2/r2 iax ), 
where nt,in is the number of galaxies in a given bin. One 
problem with this statistic, however, is that it is assum- 
ing one has measured the shear 7 instead of the reduced 
shear g, and thus overestimates the mean density when k 
is not small compared with 1. This can be corrected us- 
ing an iterative technique, similar to that used above in 
the 2-D mass reconstructions, of converting the current 



estimate of R(r) from £ c to n(r), assuming the annular 
region is at large enough radius for the subtractive ele- 
ment to be small, using this estimate of n{r) to convert g 
to 7 and recalculate £ c . As with the 2-D reconstruction, 
within 4-5 iterations the difference between successive it- 
erations of Cc are small compared to the uncertainty in 
the measurement and thus the iteration can be stopped. 
This will increase the error in the statistic due to the 
added uncertainty in the value of k used to corrected the 
reduced shear, however we have performed simulations 
of iterative aperture densitometry measurements around 
clusters with masses similar to A520 and find this ad- 
ditional error is small, usually increasing the measure- 
ment error by less than 5% compared to that calculated 
above. 

We list in Table 2 the mass measurements within a 150 
kpc radius using aperture densitometry for both statis- 
tics for each structure, in both cases using r max = 1500 
kpc and r 2 = 1150 kpc for ( c . We also show the result- 
ing mass-to-light ratios; the measurements of the clus- 
ter luminosity in these regions are slightly different than 
those given in Table 1 as we use unsmoothed luminosity 
measurements corrected for subtraction the luminosity 
density of the same annular regions used in the aperture 
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Table 2 

Aperture Densitometry Substructure Properties (r < 150 kpc) 



Substructure f Column Mass M/L (() ( c Column Mass M/L (( c ) 
(h^W 13 M Q ) (M /L zQ ) (h^W 13 M Q ) (Mq/L z q) 



PI 2.33±0.77 95±35 2.81±0.671 99±27 

P2 3.45±0.73 62±15 4.16±0.67 70±13 

P3 2.01±0.73 125±68 2.84±0.64 150±44 

P4 4.71±0.76 113±20 5.59±0.68 123±17 

P5 2.48±0.70 89±31 3.17±0.66 102±25 

P6 2.95±0.78 242±82 3.68±0.68 224±50 



Note. — M/L is calculated after subtraction of the gas mass. 



densitometry measurements. Overall the results compare 
well, within errors, with the measurements from the 2-D 
mass reconstruction given the expected increase in the 
masses in £ c due to the lack of smoothing of the aperture 
densitometry measurements, and the reduction in £ due 
to the subtraction of the mean k of the reconstruction re- 
gion. The exception is that of structure 3, for which we 
would expect the lack of smoothing to have little effect 
on the mass measurement, but instead is ~ 25% higher 
than the mass reconstruction value. These results give a 
reduced x 2 for the first 5 peaks of 1.2 for C, and 2.16 for 
£ c , which reject the constant mass-to-light ratio hypoth- 
esis at 70% and 93% confidence levels respectively. In 
both cases the main driver for the higher % 2 compared 
to that from the 2-D mass reconstruction is the lower 
mass-to-light ratio of structure 2 rather than the higher 
mass-to-light ratio of structure 3. For all three mass mea- 
surement techniques, however, a constant mass-to-light 
ratio hypothesis for all of the structures in the cluster 
cannot be excluded at a level larger than 2a. 

For completeness, we also measure the total mass of the 
cluster using the iterative aperture densitometry tech- 
nique at a radius of 700 kpc from the centroid of the 
cluster galaxy luminosity distribution, which fully con- 
tains the structures seen in Fig. [2] but is small enough to 
still have a reasonably large weak lensing S/N measure- 
ment, to measure a column mass of 5.1 ± 0.7 x 10 14 M , 
in good agreement with the values given by M07 and J12. 
Assuming that the cluster ou tside the core can still be 
mode led with a NFW profile (jN avarro. Frenk. fc W hite 
1995), using a fixed concentration c = 3.5 gives a mea- 
surement of M 2 oo = 9.1 ± 1.9 x 10 14 M©. Varying 
the concentration between 2 and 5 results in a varia- 
tion in M200 that is about half the error level of the 
weak lensing measurements. The NFW profiles have 
-M500 = 6.1 ± 1.3 x 10 14 M Q , in good agreement with 
the X-ray derived M 500 = 6.7 ± 1. x 10 14 Mq calcu- 
lated by using T x = 7.1 ± 0.7 keV dGovoni et al.ll200l 
and the M 500 - T x relation of lVikhlinin et al.1 (|2009l ). 

3.5. Strong Lensing 

One signal that could indisputably confirm the pres- 
ence of the dark peak would be a strongly lensed galaxy 
in the vicinity of region 3. A careful search for objects 
with large length-to-width ratio (giant arcs) and for ob- 
jects with similar morphologies and colors (multiply im- 
aged arclets) around region 3 in the ACS mosaic does 
not, however, reveal any obvious strongly lensed galax- 
ies. In fact, in the entire cluster system, the only two 
cases of strong lensing that we find are both being caused 




Figure 6. Shown above in a color composite from the HST ACS 
images are the only two cases of strong lensing we detect in Abell 
520. Both sets of lenses are around giant elliptical galaxies located 
in structure 2, with the strong lensing features marked by red ar- 
rows. Redshifts of the arcs are currently unknown. 

primarily by individual cluster galaxies rather than by a 
cluster core. These lenses are the two brightest elliptical 
galaxies in structure 2, and are shown in Fig. [5] The 
likelihood of having such galaxy scaled lenses is known 
to b e enhanced by the presen ce of a nearby cluster core 
(see lKneib fc N ataraian 201 1]) , however degeneracies be- 
tween the shape and strength of the galaxy's gravita- 
tional potential and that of the cluster prevent us from 
using these strong lenses to place additional constraints 
on the presence or absence of a mass peak in region 3. 
Given the cluster mass distribution's likely strong lens- 
ing cross-sections and the observed distribution, sizes, 
and redshifts of background galaxies in the field, it is not 
surprising that we do not detect any significant signs of 
strong lensing by the cluster cores. 

4. DISCUSSION 

Comparison of our results described above to those in 
J12 reveals that the weak lensing mass measurements and 
morphologies of the core structures are mostly in good 
agreement. There are, however, three major differences 
between the studies: (1) the amount of light we measure 
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Figure 7. Shown above in greyscale is what the re distribution of 
the mass reconstruction area would look like if the only signal was 
spurious ellipticity caused by the CTI of the ACS camera. The 
mass reconstruction from Fig. [2] is shown as blue contours, with 
the magenta x indicating the location of the dark peak in J12. The 
greyscale is set such that the white regions are relatively unaffected 
by the CTI, while the dark regions are the most affected. With the 
applied CTI correction, the extreme values in the CTI only mass 
reconstruction are re ~ ±0.008, which is ~ l/7th the level of the 
first contour shown in the mass reconstruction overlay. 

from cluster galaxies in location 3 is twice that measured 
in J12, (2) despite using better data, our uncertainties of 
the masses in the weak lensing measurements are 1.5 — 
1.7 times those of J12, and (3) our weak lensing column 
mass in location 3 is ~ 60 — 70% of that meaured by 
J12, depending on which measurement technique is used. 
Below we discuss possible reasons for these differences 
and potential ways to reconcile the two results. 

There are a number of differences between the inte- 
grated galaxy luminosities in our Table 1 to those in 
J12. The two largest are that ours are measured for 
the observed F814W passband (roughly a restframe R 
passband), while the J12 luminosities were measured for 
a rest-frame B passband, and our luminosities were mea- 
sured from a luminosity image that had been smoothed 
by the same amount as the mass reconstructions, while 
the J12 luminosities are measured as aperture cutouts 
without any smoothing. The smoothing decreases the 
measured luminosities in all of the structures reported in 
Table 1 , but those of structures 2 and 4 are decreased by 
a larger fraction, 30% and 35% respectively, than those 
of structures 1, 3, 5, and 6, which range from 3 — 8%. The 
column masses of the structures, however, show a similar 
increase when measured both on unsmoothed mass re- 
constructions and using the aperture densitometry tech- 
nique. As a result, while the measured mass-to-light ra- 
tio for structures 2 and 4 decrease relative to structure 
3 in the unsmoothed measurements, it is only a ~ 20% 
effect, which is not large enough to cause structure 3's 
mass-to-light ratio to be significantly different from the 
other structures. Changing from the F814W luminosi- 
ties to F435W luminosities makes the luminosity ratio of 
structure 3 to structures 2 larger and therefore moves the 
mass-to-light ratio of structure 3 closer to that of struc- 
ture 2, because all of the bright galaxies in structure 2 are 
red-sequence cluster ellipticals, while in structure 3 there 



are a number of luminous blue spiral galaxies. Looking 
at the selection of galaxies in our color selected catalogs 
reveals that all of the bright galaxies in structures 2 and 
4 are included in the catalog while several bright spiral 
galaxies were excluded by the color selection process in 
structure 3. The inclusion of these galaxies would signif- 
icantly increase the luminosity of structure 3 compared 
to structures 2 and 4. If we further reduce the size of the 
color selection region to include only cluster red-sequence 
galaxies, we do reduce the luminosity of structure 3 com- 
pared to structure 2, but only by ~ 10%. As such, the 
largest luminosity ratio between structure 2 and 3 we can 
measure is 4.1:1, compared to the 3.0:1 in Table 1, and 
can not reproduce the 5.3:1 and 6.9:1 ratios of J12 and 
M07 respectively. Thus, the majority of the difference 
between the mass-to-light ratios in the dark peak loca- 
tion between our results and those in J12 and M07 is not 
due to the weak lensing mass measurement, but comes 
from the difference in the luminosity measurements for 
galaxies at the location of structure 3. 

For the error levels in the weak lensing mass measure- 
ments in Tables 1 and 2, the two primary determining 
factors in the data are the rms shear measurement level 
and the area used to make the shear measurements. For 
the measurements in J12, the rms shear per square ar- 
cminute using their weighting scheme is 0.034 for the 
HST/WFPC2 images and 0.056 for the ground based 
CFHT image (Jee, priv comm), which are very compa- 
rable to our 0.036 from the HST/ACS images and 0.056 
from the Magellan image using our weighting scheme 
(see M3. ip . The ground-based image used in J12 is larger 
than that used here, however the effective weight for each 
galaxy in the mass measurement is 7/r oc r~ 2 , so the 
galaxies missing from our smaller field have low weight 
in the mass measurements for the 150 kpc apertures used. 
Simulations using the measured rms shear values indicate 
that doubling the size of the ground based image while 
preserving the rms shear level would decrease the size of 
the error in the mass measurements in Table 1 by only 
~ 10% of their current values. 

In addition, J12 use a much larger number density 
of background galaxies in their measurements from the 
WFPC2 images than we do from the ACS images (92 
galaxies per sq. arcminute versus our 56), and do not 
apply any color selection to their galaxies fainter than 
R = 24, but instead remove those galaxies with known 
rcdshifts that are part of the cluster and those in the red 
sequence which are brighter than R = 24. As a result, 
while they have likely excluded all of the cluster giant 
ellipticals, many cluster and foreground dwarf galaxies 
are likely still in their background galaxy catalogs as 
the majority of the galaxies in their catalogs would be at 
R > 24. With the color selection, we reject galaxies with 
an average number density of 24 per sq. arcminute (al- 
though higher towards the cluster cores than at the im- 
age edges), and using only these galaxies in a weak lens- 
ing mass reconstruction shows mostly noise with only a 
small trace of the cluster. This suggests that the majority 
of these excluded galaxies are likely cluster members or 
foreground galaxies, and the inclusion of these would de- 
crease the S /N of the weak lensing signal. The remaining 
number density of 12 galaxies per sq. arcminute differ- 
ence between the two catalogs are likely galaxies fainter 
than our minimum detection significance cut, therefore 
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having low weights in the weak lensing shear measure- 
ments and contribu ting little to the mas s reconstruction. 
By comparison, the lMiralles et all (12002ft STIS dark peak 
had a number density of 16 galaxies per sq. arcminute for 
the bright sample with the strong spurious alignment sig- 
nal, and a number density of 77 galaxies per sq. arcmin. 
giving the weaker alignment signal. Overall, the inclu- 
sion of faint cluster and foreground galaxies will increase 
the value of £ C rit needed to convert k measurements into 
a surface mass, which is consistent with the J 12 adopted 
E crit = 4.1 x 10 9 M Q /kpc 2 versus our E crlt = 3.6 x 10 9 
M Q /kpc 2 . This higher value of S C rit gives a higher error 
in the surface mass for a given error in n. Simulating the 
reported error levels for the shear measurements in J12 
with the reported £ cr i t suggests that they should have 
error levels on the order of ~ 7 x 10 12 M Q rather than 
the ~4x 10 12 M Q level that they report, and thus their 
significance levels are less than 60% of what they state. 

The most surprising difference between our results and 
those of J 12, however, is the similarity of the weak lens- 
ing measurements for all parts of the cluster except that 
in the vicinity of region 3. When an error is made in the 
analysis of a weak lensing data set, the resulting mass 
reconstruction normally differs from the true mass re- 
construction across the entire field, and not a localized 
difference in one small portion of the mass reconstruc- 
tion, especially if the difference in that one area is large 
compared to the noise level in the reconstructions. Both 
J12 and M07 overstated the significance of their detec- 
tions of this structure because they compared their mea- 
surements to a mass level in that region instead of to a 
model of a constant mass-to-light ratio across the cluster. 
We argue that our mass measurement of 2.3 x 10 13 M Q 
for this region is a better baseline for comparison as it 
has a similar mass-to-light ratio as the other structures, 
and this would give a significance of the mass overden- 
sity in this region of 1.9a for the M07 measurement and 
4ct for the J12 measurement, using their quoted errors, 
or around 2.3cr using our estimated error level for their 
measurement discussed above. The significance of the 
difference in the mass measurements for this region be- 
tween our result and the J12 result is larger than the 2. So- 
level, however, as one of the primary sources of noise, the 
intrinsic ellipticity of the background galaxies used in the 
weak lensing shear measurement, is largely in common 
between the two measurements. Any deviation in shear 
measurements between the two catalogs would therefore 
either have to be an error in one of the measurement 
techniques, or a much more significant inherent align- 
ment in the shapes of the galaxies used in J12 that are 
not included in our background galaxy catalogs. At this 
level of difference, there are plenty of potential sources of 
error in the weak lensing measurements to explain how 
two such different reconstructions can be drawn from the 
same underlying data. 

One potential source of error in a weak lensing mass 
reconstruction is that arising from an incorrect PSF 
smearing correction. This is cited in M07 to explain 
the difference in their reconstruction of the Subaru 
im age compared wi t h an analysis in a early version 
of lOkabe fc Umetsul (|2008f ). which Okabe and Umetsu 
agreed with in the final version of their paper. It is un- 
likely an error in PSF correction is a factor for this case. 



Both groups have tested their methodology on simula- 
tions and recover shears to better than a 1% accuracy 
with known PSFs, well below the level that would create 
or remove such a large structure in the mass reconstruc- 
tion. Also, both groups have multiple image sets from 
different telescopes on which they have measured similar 
looking structures, making the likelihood of the presence 
or absence of a given structure being caused by a poorly 
measured PSF in that region small. Both J12 and M07 
find a mass structure in region 3 using two ground based 
image sets from different telescopes and, for J12, a HST 
WFPC2 mosaic. We show in Fig. [5] separate mass recon- 
structions for the Magellan ground based image and the 
HST ACS mosaic, neither of which exhibit a structure 
in region 3. Because the different imaging sets have very 
different intrinsic PSF ellipticity distributions, it would 
be highly unlikely to make mistakes in each shear mea- 
surement such that they combine to give the same false 
weak lensing mass structure. 

Another potential source of error is the treatment of 
the charge transfer inefficiency (CTI) in the HST im- 
ages. Both WFPC2 and ACS suffer from CTI, which is 
caused primarily by electron traps being created in the 
silicon detectors by high energy cosmic rays. As these 
traps effectively shuffle charge to trails behind a bright 
object, they impart a fake, correlated ellipticity in the 
stars and galaxies that can mimic a weak lensing signal. 
The WFPC2 CTI shape was measured by J12 and while 
they do find a potential residual in the location of the 
dark peak location, the residual is more than an order 
of magnitude below what would be needed to cause the 
observed structure. For the ACS images, we are con- 
fident that CTI is not playing a significant role in the 
galaxy and PSF shape measurement process for two rea- 
sons. The first is that we created two different sets of 
images using different CTI correction mechanisms (see 
H2.1[) and measure nearly identical mass reconstructions. 
The second is that the induced ellipticity pattern from 
the ACS CTI when used in the 2-D mass reconstruction 
technique would create a four quadrant positive k signal 
as well as significant k peaks both north and south of the 
HST imaging area which is not seen in our A520 mass 
reconstruction (see Fig. [7]). If both of the CTI correction 
techniques were ineffective in removing the induced shear 
signal, correction for this would not significantly change 
the mass measurement for region 3 as our tiling strategy 
for ACS placed it in a region with no CTI. However, it 
would reduce the masses of all of the other structures by 
25 — 40%, placing these masses in disagreement with that 
measured by J12 and that which we measure from only 
the Magellan data. This would also shift the centroids 
of the weak lensing peaks for structures 1, 2, 4, and 5 
away from the black bars in Fig. [71 which would cause 
significant offsets between the mass and galaxy locations 
in these structures. The reduction in the CTI effect by 
a factor of 20 from the correction techniques (§ I2.1j) will 
reduce the CTI caused k features by also roughly a factor 
of 20, which would have a maximum change of ~ l/7th 
the value of the first contour of the overlays and change 
the measured masses in Table 1 by order 1 — 2% for all 
but structure 3, which would still be unaffected by the 
CTI. 

A third potential source of error is a sudden change in 
the mean shear at the boundary of the ACS images due 
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to differences in the redshift distributions of the back- 
ground galaxies used in the ACS images compared to 
those in the Magellan image. This would cause a sud- 
den change in the mean £ cr j t that, as the mass recon- 
struction algorithm assumes S C rit is constant across the 
field, would be interpreted as a feature caused by the 
surface mass distribution of the cluster. In simulations, 
such a mismatch usually causes a plateau of surface den- 
sity with a size and shape similar to the ACS mosaic 
region, but can occasionally due to the coupling of ran- 
dom noise in the shear measurements and this systematic 
change at the boundary alter the relative masses and lo- 
cations of the structures in the cluster core seen in the 
mass reconstructions. We already attempted to correct 
for this effect by scaling the shear measurements of the 
ground based image to match the depth of the ACS im- 
age prior to combining the catalogs. To test that this is 
not causing a problem, we have recombined the space- 
based shear catalog with both an unsealed ground-based 
catalog and one with twice the correction we originally 
used, and with neither of the two resulting catalogs does 
the mass reconstruction have any structure at the dark 
peak location (3). We note that this type of rescaling 
of the ground-based shear measurements by a constant 
fraction is not strictly correct, as we are really measuring 
the reduced shear and both k and 7 will change with a 
variation in E cr j t . As a final test, we used the original 
mass reconstruction to separate the 7 and k components 
of the reduced shear in the ground based image cata- 
log, scaled each separately, and recombined them to get 
a scaled reduced shear catalog. Combining this catalog 
with the ACS catalog also made no discernible change in 
the masses or centroids of the structures observed in the 
resulting mass reconstruction. 

A more likely cause of the difference between the two 
reconstructions is not an error in the weak lensing mea- 
surements at all, it is simply a difference in which galaxies 
are selected to have their shapes measured and included 
in the shear field. Because the intrinsic galaxy shapes 
are the largest random error for reasonably bright galax- 
ies, different sets of galaxies will change the noise field 
in the reconstructions. It is this source of noise that is 
reproduced in the bootstrap resampling methodology of 
t j3.3l As can be seen in Fig. [31 we do see a structure 
near location 3 in about 2% these bootstrap resampled 
reconstructions, which is consistent with our excluding 
the mass measured by J12 in this region at a 2.3ct level. 
The additional cluster and foreground galaxies in the J12 
catalogs discussed above could be a source for the addi- 
tional shear needed to produce a peak in location 3, ei- 
ther through intrinsic or purely random alignments of the 
intrinsic shapes of the galaxies s uch as was seen in th e 
"dark lens" of the STIS fields in iMiralles et~aTl (pOOl . 
To know for certain the cause, however, will require a di- 
rect galaxy by galaxy comparison of the two catalogs to 
determine if the difference is within the shear measure- 
ments for galaxies in common between the two catalogs 
or from the extra galaxies used in the J12 catalog. 

There is one feature in common in the mass recon- 
structions which is seemingly at odds with a CDM model 
of the merger, however. East of structure 5 are several 
bright elliptical galaxies which have sp e ctroscopic red- 
shifts dYee. Ellingson. fc Carlberd Il996t iCarlberg et~aTl 
Il996t iProust et al.ll2000|) that identify them as cluster 



members, but neither our nor the J 12 mass reconstruc- 
tion show any appreciable mass overdensity around these 
galaxies despite them being several times more luminous 
that those of structure 5. The bootstrap resampled cat- 
alogs do have peak 5 moving far enough eastward to be 
associated with these galaxies, but only ~ 10% of the 
time. If real, this lack of mass could be an example of 
ejection of bright galaxies from cores during a merger 
event. However, as these galaxies are located at the edge 
of both the ACS and WFPC2 mosaics, it is possible that 
a boundary effect problem, as discussed above, could be 
causing a drop in surface density at this location. In the 
Magellan only reconstruction seen in Fig. [SJ we do find 
structure 5 to have an extended mass tail across this lo- 
cation, although the centroid is still more consistent with 
the fainter ellipticals of location 5. 

5. CONCLUSIONS 

Using a new multi-color HST ACS mosaic and previ- 
ously unpublished Magellan image set, we performed a 
weak lensing mass reconstruction on the merging cluster 
A520. The mass structures in the reconstruction show 
excellent agreement with the distribution of light from 
cluster galaxies after subtraction of the mass of the intra- 
cluster X-ray plasma. While the masses we measure for 
the cluster overall and all of the cluster substructures 
containing galaxies are in good agreement with previ- 
ous weak lensing measurements in J12, we do not detect 
the mass overdensity spatially coincident with the X-ray 
plasma cloud that was found in both M07 and J12. We 
measure a total mass in this region consistent with a con- 
stant mass-to-light ratio across the cluster, and exclude 
the additional mass in the central region at a ~ 98% 
confidence level. Using an aperture densitometry mea- 
surement instead of the mass reconstruction results in a 
slightly higher mass for the dark peak region, it still ex- 
cludes the mass from J12 at a 93% confidence level. The 
mass measurements from the 2-D mass reconstruction 
are consistent with a constant mass-to-light model, while 
the mass aperture show marginal evidence for a depar- 
ture from constant mass-to-light ratio, primarily caused 
by a lower mass-to-light ratio measurement than aver- 
age in structure 2. We also find that the significances 
for the dark peak structure were overstated in both M07 
and J12 as they calculated the significance by compar- 
ing their measured mass to a mass of in the center of 
the cluster, and their significances of detection are < 2a 
and ~ 2.3cr respectively when measured compared to our 
constant mass-to-light ratio model. 

We have considered several potential causes for the dis- 
crepancy in the central region between the various mass 
reconstructions while still having good agreement in the 
rest of the cluster core. We suggest that the most likely 
explanation is an inherent alignment in the galaxies that 
were included in the J12 shear measurements but ex- 
cluded from ours. Regarding the discrepent mass-to- light 
ratio, we find that both M07 and J12 have significantly 
lower cluster luminosity measurements in the region of 
the "dark peak" than our measurements, and this differ- 
ence in cluster luminosity is responsible for the majority 
of the difference in the mass-to-light ratios for the central 
region between the studies. 

We identify one structure on the eastern edge of the 
HST image which has bright elliptical galaxies that are 
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known to be part of the cluster for which neither we nor 
J12 obtain a significant amount of mass. However, we 
do detect mass in this region in a mass reconstruction 
shape measurements from the Magellan image. We are 
uncertain if the lack of mass in this region is an aspect 
of the merger or an edge effect from the HST mosaic. 

The overall mass structure that we measure for A520 
is in good agreement with a constant mass-to-light ratio, 
and therefore with collisionless cold dark matter — sim- 
ilar to the conclusions drawn from all other well-studied 
merging clusters. Deriving a quantitative upper limit 
on the dark matter self-interaction cross-section from 
A520 will require additional kinematic information and 
detailed modeling of this merging system. 
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